gD = ggplotGrob(sigma)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB,  size = "last"),
gshape
), size = "first"
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB,  size = "last"),
gShape
), size = "first"
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
gShape, ncol = 1
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
gShape, ncol = 1
)
dev.off()
setwd(dir)
source("mod_hierarchy_shape.R")
getwd()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
gShape, ncol = 1
)
dev.off()
setwd(dir)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
gShape, ncol = 1
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange( gA, gB, gShape)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange( gA, gB, gShape, ncol = 2)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
gShape, ncol = 1
)
dev.off()
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_hierarchy_shape.R')
firm.shape = ggplot() +
geom_bar(data = mod, aes(x = level, y = median),  fill = "dodgerblue4", alpha = 0.6, stat = "identity",  position = "dodge", width = 0.94) +
geom_errorbar(data = mod, aes(x = level, ymin = lwr, ymax = upr), width = 0.2) +
scale_x_continuous("Hierarchical Rank", breaks = seq(1, 12, 1)) +
scale_y_continuous("Percent of Employment", breaks = seq(-100, 100, 10), labels = abs(seq(-100, 100, 10)) ) +
coord_flip(xlim = c(0.5, 11)) +
ggtitle("C.  Shape of a Generic Firm in the Hierarchy Model") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(color = "black"),
text=element_text(size=text.size, family="Times"),
axis.title=element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt"))) +
annotate("text", label = "Variation in Rank Size", x = 2.5, y = 38, family = "Times", size = 2.8) +
geom_segment(aes(x = 1.95, y = 38, xend = 1.2, yend = 38),  arrow=arrow(type="closed", length=unit(0.1,"cm")), color = "grey50")
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
gA = ggplotGrob(a)
gB = ggplotGrob(b)
gShape = ggplotGrob(firm.shape)
gC = ggplotGrob(gini.hist)
gD = ggplotGrob(sigma)
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB),
cbind(gShape, blank, size = "first", widths = c(2,1))
)
)
dev.off()
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB),
cbind(gShape, blank, size = "first")
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first", widths = c(2,1))
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first")
)
)
dev.off()
gA = ggplotGrob(a)
gB = ggplotGrob(b)
gShape = ggplotGrob(firm.shape)
gC = ggplotGrob(gini.hist)
gD = ggplotGrob(sigma)
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first")
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first"),
size = "first"
)
)
dev.off()
?widthDetails
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first", ),
size = "first"
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first", widths = c(3,1) ),
size = "first"
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first"),
size = "first"
)
)
dev.off()
library(egg)
png("span_parameters.png", width = 7, height = 6, res = 600, units = "in")
ggarrange(gA, gB, widths = c(1,2))
dev.off()
png("span_parameters.png", width = 7, height = 6, res = 600, units = "in")
ggarrange(a, b, widths = c(1,2))
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
ggarrange(a, b, widths = c(1,2))
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gShape, blank, size = "first"),
size = "first"
)
)
dev.off()
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
arrangeGrob(gShape, blank, ncol = 2, widths = c(3,1)),
ncol = 1
)
dev.off()
firm.shape = ggplot() +
geom_bar(data = mod, aes(x = level, y = median),  fill = "dodgerblue4", alpha = 0.6, stat = "identity",  position = "dodge", width = 0.94) +
geom_errorbar(data = mod, aes(x = level, ymin = lwr, ymax = upr), width = 0.2) +
scale_x_continuous("Hierarchical Rank", breaks = seq(1, 12, 1)) +
scale_y_continuous("Percent of Employment", breaks = seq(-100, 100, 10), labels = abs(seq(-100, 100, 10)) ) +
coord_flip(xlim = c(0.5, 11)) +
ggtitle("C.  Shape of a Generic Firm in the Hierarchy Model") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(color = "black"),
text=element_text(size=text.size, family="Times"),
axis.title=element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt"))) +
annotate("text", label = "Variation in\nRank Size", x = 3, y = 38, family = "Times", size = 2.8) +
geom_segment(aes(x = 1.95, y = 38, xend = 1.2, yend = 38),  arrow=arrow(type="closed", length=unit(0.1,"cm")), color = "grey50")
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
arrangeGrob(gShape, blank, ncol = 2, widths = c(3,1)),
ncol = 1
)
dev.off()
firm.shape = ggplot() +
geom_bar(data = mod, aes(x = level, y = median),  fill = "dodgerblue4", alpha = 0.6, stat = "identity",  position = "dodge", width = 0.94) +
geom_errorbar(data = mod, aes(x = level, ymin = lwr, ymax = upr), width = 0.2) +
scale_x_continuous("Hierarchical Rank", breaks = seq(1, 12, 1)) +
scale_y_continuous("Percent of Employment", breaks = seq(-100, 100, 10), labels = abs(seq(-100, 100, 10)) ) +
coord_flip(xlim = c(0.5, 11)) +
ggtitle("C.  Shape of a Generic Firm in the Hierarchy Model") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(color = "black"),
text=element_text(size=text.size, family="Times"),
axis.title=element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt"))) +
annotate("text", label = "Variation in\nRank Size", x = 3, y = 38, family = "Times", size = 2.8) +
geom_segment(aes(x = 1.95, y = 38, xend = 1.2, yend = 38),  arrow=arrow(type="closed", length=unit(0.1,"cm")), color = "grey50")
gA = ggplotGrob(a)
gB = ggplotGrob(b)
gShape = ggplotGrob(firm.shape)
gC = ggplotGrob(gini.hist)
gD = ggplotGrob(sigma)
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
arrangeGrob(gShape, blank, ncol = 2, widths = c(3,1)),
ncol = 1
)
dev.off()
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
b = ggplot() +
geom_density(data = model, aes(b), fill = "dodgerblue3", alpha = 0.5, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1, 0.025)) +
labs(x = expression(italic(b)), y = "Density") +
ggtitle(bquote(bold(B.~~Model~Parameter~italic(b)))) +
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(),
text=element_text(size=text.size, family="Times"),
axis.title.x = element_text(size=rel(0.97)),
axis.title.y = element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt")))
b = ggplot() +
geom_density(data = model, aes(b), fill = "dodgerblue3", alpha = 0.5, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1, 0.025)) +
labs(x = expression(italic(b)), y = "Density") +
ggtitle(bquote(bold(B.~~Model~Parameter~italic(b)))) +
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(),
text=element_text(size=text.size, family="Times"),
axis.title.x = element_text(size=rel(0.97)),
axis.title.y = element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt")))
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_hierarchy_shape.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
setwd(dir)
gA = ggplotGrob(a)
gB = ggplotGrob(b)
gShape = ggplotGrob(firm.shape)
gC = ggplotGrob(gini.hist)
gD = ggplotGrob(sigma)
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
arrangeGrob(gShape, blank, ncol = 2, widths = c(3,1)),
ncol = 1
)
dev.off()
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
a = ggplot() +
geom_density(data = model, aes(a), fill = "red", alpha = 0.5, adjust = 2) +
scale_y_continuous("Density", breaks = c(0:4), labels = c("~0", " 1", " 2", " 3", " 4")) +
labs(x = expression(italic(a))) +
ggtitle("A.  Variation in Parameter a") +
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(),
text=element_text(size=text.size, family="Times"),
axis.title.x = element_text(size=rel(0.97)),
axis.title.y = element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt")))
a = ggplot() +
geom_density(data = model, aes(a), fill = "red", alpha = 0.5, adjust = 2) +
scale_y_continuous("Density", breaks = c(0:4), labels = c(" 0", " 1", " 2", " 3", " 4")) +
labs(x = expression(italic(a))) +
ggtitle("A.  Variation in Parameter a") +
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(),
text=element_text(size=text.size, family="Times"),
axis.title.x = element_text(size=rel(0.97)),
axis.title.y = element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt")))
library(ggplot2)
library(readr)
library(gridExtra)
library(grid)
library(gtable)
library(data.table)
library(hmod)
library(snow)
library(doSNOW)
library(foreach)
library(tcltk)
library(here)
library(boot)
n.iteration = 10000
dir = here()
wd = gsub("Figures", "Hierarchy Model/data", dir)
setwd(wd)
s = fread("span_empirical.txt")
g = fread("within_level_gini_empirical.txt")
cl = makeCluster(8, type="SOCK")
registerDoSNOW(cl)
clusterSetupRNG (cl, type = "RNGstream")
# Progress Bar
beginTime = proc.time()
pb = txtProgressBar(max = n.iteration, style = 3)
progress = function(n) setTxtProgressBar(pb, n)
opts = list(progress=progress)
# 4. Run model for each set of parameters  ----------------------------------------------------
model =  foreach ( iteration = 1:n.iteration, .options.snow=opts, .combine = rbind, .packages=c('hmod')) %dopar%  {
span = boot_span(s$V1, s$V2)
sigma = boot_sigma(g$V1)
c(span, sigma)
}
stopCluster(cl)
model = data.table(model)
names(model) = c("a", "b", "sigma")
# plots
##########################################
text.size = 10
a = ggplot() +
geom_density(data = model, aes(a), fill = "red", alpha = 0.5, adjust = 2) +
scale_y_continuous("Density", breaks = c(0:4), labels = c(" 0", " 1", " 2", " 3", " 4")) +
labs(x = expression(italic(a))) +
ggtitle("A.  Variation in Parameter a") +
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
axis.line = element_line(),
text=element_text(size=text.size, family="Times"),
axis.title.x = element_text(size=rel(0.97)),
axis.title.y = element_text(size=rel(0.97)),
plot.title = element_text(size=rel(1), face="bold", hjust = 0.5),
axis.ticks.length = unit(-0.7, "mm"),
axis.text.x = element_text(margin=margin(6,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,2,3,"pt")))
setwd(dir)
gA = ggplotGrob(a)
gB = ggplotGrob(b)
gShape = ggplotGrob(firm.shape)
gC = ggplotGrob(gini.hist)
gD = ggplotGrob(sigma)
blank <- ggplot() + theme_void()
blank = ggplotGrob(blank)
pdf("span_parameters.pdf", width = 7, height = 6)
grid.arrange(
arrangeGrob(gA, gB, ncol = 2),
arrangeGrob(gShape, blank, ncol = 2, widths = c(3,1)),
ncol = 1
)
dev.off()
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/case_parameters.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/cap_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/q_q_plot.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/cap_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_table.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/cap_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/q_q_plot.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_table.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/labor_mod_properties.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
plot(unlist(h_power[,1]))
plot(unlist(h_power[1,]))
plot(unlist(h_power[1,]), log = "y")
plot(unlist(h_power[3,]), log = "y")
plot(unlist(h_power[5,]), log = "y")
plot(unlist(h_power[20,]), log = "y")
plot(unlist(h_power[50,]), log = "y")
plot(unlist(h_power[30,]), log = "y")
plot(unlist(h_power[80,]), log = "y")
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
plot(unlist(h_power[80,]), log = "y")
plot(unlist(h_power[30,]), log = "y")
plot(unlist(h_power[80,]), log = "y")
plot(unlist(h_power[40,]), log = "y")
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
plot(unlist(h_power[80,]), log = "y")
plot(unlist(h_power[30,]), log = "y")
plot(unlist(h_power[20,]), log = "y")
plot(unlist(h_power[10,]), log = "y")
plot(unlist(h_power[15,]), log = "y")
plot(unlist(h_power[16,]), log = "y")
plot(unlist(h_power[65,]), log = "y")
plot(unlist(h_power[69,]), log = "y")
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
plot(unlist(h_power[69,]), log = "y")
plot(unlist(h_power[65,]), log = "y")
plot(unlist(h_power[16,]), log = "y")
plot(unlist(h_power[65,]), log = "y")
plot(unlist(h_power[66,]), log = "y")
plot(unlist(h_power[69,]), log = "y")
plot(unlist(h_power[90,]), log = "y")
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/k_frac_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/h_power_percentile.R')
source('~/Dropbox/cloud_work/papers in progress/rich_different_2/Supplementary Materials/Figures/mod_properties.R')
install.packages("ks")
source('~/Dropbox/cloud_work/papers in progress/rich_different_3/Supplementary Materials/Figures/beta_regress.R')
